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Abstract 

A general problem in optimal control consists of finding a terminal 
reward that makes the value function independent of the horizon. Such 
a terminal reward can be interpreted as a max-plus eigenvector of the 
associated Lax-Oleinik semigroup. We give a representation formula for 
all these eigenvectors, which applies to optimal control problems in which 
the state space is non compact. This representation involves an abstract 
boundary of the state space, which extends the boundary of metric spaces 
defined in terms of Busemann functions (the horoboundary) . Extremal 
generators of the eigenspace correspond to certain boundary points, which 
are the limit of almost-geodesics. We illustrate our results in the case of 
a linear quadratic problem. 



1 Introduction 

A basic optimal control problem with finite horizon consists of maximising a 
reward of the form 

/ L(x(s),u(s))ds + 0(x(t)) , (la) 

^0 

over the trajectories x(-) and u(-) satisfying 

x(s)-/(x(s),u(s)) , (lb) 
u(s) e U, x(s) £ X for all < s < t, (Ic) 
and x(0) = X . (Id) 

We are given the state space X C M", the set of admissible controls U C M.^, 
the dynamics f : X x U ^ R", the initial condition x E R", the horizon t, the 
Lagrangian L : X x U M, and the terminal reward : X ^ R U {— oo}. Of 
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course, the data of the problem must satisfy certain technical assumptions that 
we skip for the moment, and the trajectories x(-) and u(-) must have a certain 
regularity: u is required to be measurable and x is required to be absolutely 
continuous. We are interested in the maximal value of the reward, which we 
denote by ^'(x), and in the optimal solutions x(-) and u(-). 

One may wish in particular to determine a final reward 4> making the value 
function : X — > M U {— cx)}, x i-^ v*{x), independent of the horizon i in a 
reasonable sense. Since u* is expected to grow linearly with the horizon i, we 
look for final rewards such that 



where A is a constant that can be interpreted as a mean reward per time unit. 
(For any function and constant /i, we denote by /i + the function sending x 
to fj. + (j){x).) 

In this paper, we address the problem of finding all the solutions (f) oi ^ 
corresponding to a given constant A. 

The economic motivation of this classical problem is presented in Pj: in 
finite horizon optimal control, the optimal strategy may change dramatically as 
the current time gets closer to the horizon, which is often undesirable. To give 
a simple illustration, imagine a company having a concession to run a network 
belonging to a city during a given period. The company might stop repairing the 
network as time approaches the horizon. Thus, the city might wish to choose a 
terminal penalty (p preventing the company's strategy from being distorted by 
the horizon. A way of addressing this issue is to solve (0). 

To study this problem, it is convenient to introduce the Lax-Oleinik semi- 
group^ {S*)t>Q, which consists of the operators 5* sending the terminal reward 
to the value function w* . Observe that the operator S** is max-plus linear, mean- 
ing that it preserves suprema and commutes with the addition of a constant. 
Equation 10) can be rewritten as 



We recognise this as a max-plus analogue of the spectral equation. We say that 
(j) is an eigenfunction of the semigroup and that A is the associated eigenvalue. 

The discrete-time analogue of problem JJJ has received a considerable amount 
of attention 01 E E]. Problem (jSl itself has been studied by Maslov, 
Kolokoltsov, and other members of the "idempotcnt analysis" school |3IH]. More 
recently, it has appeared in the setting of Fathi's "weak KAM" theory [Hl llOlITT] . 
Results of idempotcnt analysis and of weak KAM theory yield a representation 
of all the solutions of Q, under certain assumptions. For instance, the ap- 
proach of II is in the setting where X is a Riemannian manifold and the 
Lagrangian has smoothness and strict-convexity properties. The representation 
results of [3 Section 3.2], on the other hand, apply to other special cases. How- 
ever, the discrete-time theory 6 suggests that general results should hold in the 
continuous-time setting. Such a generality is required in many applications. For 
instance, the Lagrangian might not be smooth or the trajectory x may be the 
sum of an absolutely continuous part and a singular part modelling "jumps" . 

In , some continuous-time results were obtained in special cases by apply- 
ing the discrete-time results. However, we obtain more general results when we 
redo the theory from the start in continuous-time, which we do here. Our main 



V* = Xt + (p for alH > , 



(2) 



S*(j) ^Xt + (j) for alH > . 



(3) 
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results, Theorems l3. 15l and l3. 1 II show that the solutions of Q that are extremal 
in the max-plus sense can be identified with limits of "almost-geodesics" , and 
that any solution is a max-plus linear combination of extremal solutions. Our 
study relies on a compactification of the state space X, which is similar to the 
compactification of metric spaces by horofunctions (generalised Busemann func- 
tions) jl2[ [T5| and to the Martin compactification in potential theory |14[ [T5] . 

The interest of the notion of almost-geodesic is that it leads to the most 
general results since it allows us to avoid making hypotheses on the Lagrangian 
or the space X that imply the existence of geodesies. 

In Sectional we illustrate these results by discussing a simple optimal control 
problem. We take the state space X to be R", the set of controls U to be M", 
and /(x, u) := u. The Lagrangian is 

L(x,u):^~\x\^ -\u\^ , (4) 

where | • | denotes the Euclidean norm. 

This optimal control problem is linear quadratic. In the case where A = 0, it 
is well known that there are quadratic solutions of Q , which can be computed 
by solving the algebraic Riccati equation. These solutions are = ±|a;p. The 
solution in which one is traditionally interested is the "stable" one (f){x) = — \x\'^. 
It yields the feedback control u{t) = V0(x(i)) = — 2x(t), which pushes the state 
towards the origin. The other quadratic solution yields the opposite control, 
sending the state off to infinity, radially. We show however that there are other 
solutions, which are non-quadratic. Amongst them, the extremal solutions, like 
the unstable quadratic one, send the state off to infinity, each in a different 
direction. 

When A > 0, there is nothing like the stable quadratic solution, but the 
non-quadratic solutions persist. 

2 The max-plus Martin space of a max-plus lin- 
ear semigroup 

Recall that the max-plus semiring Rmax is the set RU {— oo}, equipped with the 
addition operation (a, b) t-^ max(a, b) and multiplication operation (a, b) i— s- a+b. 
The completed max-plus semiring Rmax is obtained by adjoining to Rmax a 
-|-oo element, with the convention that — oo is absorbing for (a, b) i— > a + b. 
Although in optimal control applications, the state space X is typically a subset 
of R", it is convenient to assume only that X is an arbitrary set. To any map 
B : X X X Rmax, (a;,?/) i-^ B{x,y) is associated a max-plus linear self-map 
from Rmax to itself, g i— s- Eg, defined by: 

Bg{x) = sup {B{x, y) + g{y)) . 
yex 

We say that (x, y) i—)- B(x, y) is the kernel of the map g i—s- Bg. It is uniquely 
defined. It is known that a large class of max-plus linear maps can be represented 
in this form, see |16j . 



Definition 2.1. We say that a family of operators (A*)t>o from Rmax to itself 
is a max-plus linear semigroup with kernel if every A* can be written as 

A'g{x)^sn^{A\x,y)+g{y)) , 
yex 
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for some function {x,y) i-^ A^(x,y) G Mmax, and if 

A° = 1 and A'+^ = A'A'^ , for all s, t > , 

where concatenation denotes the composition of operators, and / denotes the 
identity operator. 

The semigroup property, = A^A*, is equivalent to 

A^+\x,y) = snp{A^x,z)+A\z,y)) . (5) 

Example 2.2. The Lax-Oleinik semigroup is an example of a max-plus linear 
semigroup with kernel: S^{x,y) is the supremum of pa|l over all trajectories x 
and u satisfying ljlb|l - l)ldp and x(i) = y. In this case the semigroup property 
follows from the dynamic programming principle. 

If (A*)t>o is a max-plus linear semigroup with kernel, then for all A G M, 

Ai: g^-\t + A'g 

also defines a max-plus linear semigroup, with kernel ^^(a;, y) — ~Xt + A'^{x, y). 
Therefore, when looking for the solutions of it is legitimate to assume 
that A = 0. The analogy with potential theory justifies the following definition. 

Definition 2.3. We say that a function (f) : X ^ Kmax is harmonic with respect 
to the semigroup (>l*)t>o if 

A*0 = (j) for alH > . 

We say that 4> is super-harmonic if 

< (j) for alH > . 

We denote by and y the sets of harmonic and super-harmonic functions, 
respectively. 

Observe that (p is required not to take the +oo value. We define 

A* := sup A* . 
t>o 

(The sup is for the pointwise ordering of operators.) To study harmonic func- 
tions, we need to exclude degenerate cases. 

Assumption 2.4. We assume that A*{x, y) is finite for all x and y in X . 

Observe that A* A* ~ A* , and so 

A*{x,z) + A*{z,y) < A*(x,y) for aU x, z G X. (6) 

This is an analogue of the triangular inequality. By taking x — y — z, we deduce 
that A*{x,x) < 0. Since A* > I, we have A*{x,x) > 0, and so A*{x,x) = 0, for 
all X £ X . We fix an arbitrary point b £ X, which we call the hasepoint. 
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Definition 2.5. The (max-plus) Martin kernel of the semigroup {A*)t>o with 
respect to the basepoint b is defined by: 

K{x,y)=A*ix,y)-A*{b,y) . 

The (max-plus) Martin space of the semigroup (A*)(>o is the closure in the 
topology of pointwise convergence of the set := {K{-,y^ \ y e X\ C M'^. 
The (max-plus) Martin boundary is £^ ~ ^ \ . 

From linj, we deduce that A*{b,x) + A*(x,y) < A*{b,y), hence 

Kix,y) < -A*{b,x) , for all a;, 2/ e X . 

Similarly, 

K{x, y) > A*{x, b) , for al\x,y e X . 

Thus C ]^^gj^[A*(a;, 6), — A*(6, x)], and the latter space is compact, by Ty- 
chonoff's theorem. Hence, the closure ^ of is compact. 

Remark 2.6. In many applications, the state space X is naturally equipped 
with a metric d. For instance, when X C K", we can take for d the metric 
induced by any norm on R". Then, the Martin compactification could be defined 
alternatively with respect to the topology of uniform convergence on compact 
sets, or on bounded sets. This often leads to the same Martin compactification. 
Indeed, the triangular inequality yields 

A* (x, y) < K{x, z) - K{y, z) < -A* (y, x) . (7) 

It is often the case that |A*(x, y)\ < Cd{x, y) for some constant C > 0. Then, {T)) 
shows that the set J(f is equicontinuous, and so, by Ascoli's theorem, its closure 
in the topology of pointwise convergence is the same as its closure in the topology 
of uniform convergence on compact sets. 

Any element of is super-harmonic, because A* A* < A* . Since a pointwise 
limit of super-harmonic functions is super-harmonic, we deduce the following. 

Proposition 2.7. Each element of .y^ is super-harmonic. 

As in the discrete time case 0, we have the following. 

Proposition 2.8. We have K{-, x) — K{-, y) if and only if 

A*{x,y)+A*{y,x)^0 . (8) 

In the special case of the Lax-Oleinik semigroup, JHJ means that we can find 
trajectories starting from x, passing to y, and returning to x, with a reward 
arbitrarily close to 0. 

Define the equivalence relation on X such that x ^ y ii ^ holds. Propo- 
sition allows us to identify J(f with the quotient X X/ ^. In particular, 
when the equivalence classes of ^ are singletons, the map x t-^ K{-,x) may 
be thought of as an "embedding" of X in .4^, and the elements of the Martin 
boundary — ^ \ Ji^ may be thought of as "boundary points" of X. 

For all functions X ^ Kmax and for all r\ G we set: 

fJ-iiv) limsup {A*{b,x) + ^{x)) 

:= inf sup {A*{b,x) + ^{x)) , 
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where the infimum is taken over all neighbourhoods of 77 in ^ . This map 
will play the role of the spectral measure in the Martin representation theorem. 
As in the discrete time case, we have the following. 

Proposition 2.9. Let u be a super-harmonic function. Then, ^u{K{-^x)) — 
A*(h,x) -\- u{x) for each x E X. 

This result means that the map X — s- Mmax, 2; i-^ A*{h, x) + u{x) induces a 
map on X, or equivalently on J>^, which is upper-semicontinuous in the topology 
on J^f and that /i„ extends this map. 

It is convenient to renormalizc the operator A* by means of the following 
transformation. We define, for all x,y G X, 

A\x,y) = A*{b,x)+A*{x,y) - A*{b,y) . 

The triangular inequality implies that 

A^{x,y) < 0, for al\x,yeX . 

Proposition 12 . 91 allows us to extend the kernel A^ to a kernel H on x ^ 
as follows: 

Indeed, 

H{K{;x),K{;y))^A\x,y) . 

Definition 2.10. We define the minimal boundary ^™ to be the set of points 
5 of ^ that are (max-plus) harmonic and satisfy C) = 0- 

3 Main results 

In this section, we state our main results. Many of the proofs are similar to those 
of where discrete time semigroups are considered. However, the notion of 
almost-geodesic appearing there must be adapted to continuous time. 

We use the term path for an arbitrary map 7 from a closed interval of M+ 
to X. We set, for alHg < < • • • < in the interval of definition of 7, 

n-l 

J(to, . . . , i„; 7) := 5Z l{U+i)). 

4=0 

Definition 3.1. An almost-geodesic with parameter e is a path 7 : K+ —^ X 
such that, for all n G N and < ii < i2 < • • • < ^n, we have 

/(ii,...,i„;7)>^*(7(ii),7(in))-e- 

Definition 3.2. A path 7 : M+ ^ X is said to be e- almost-optimal with respect 
to a super harmonic function h if for all n G N and = < ti < t2 < • • • < tn, 
we have 

< e + I{to, . . . , t„; 7) + h{j{t^)). 

Lemma 3.3. Let 7 : IR+ X be an almost-optimal path with respect to a 
superharmonic function. Then 7 is an almost-geodesic. 
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Definition 3.4. The reward of a path 7 : [sq, si] ^ X is 
^(7) := , inf /(to, ■ • ■,tn;i), 

where the infimum is taken over all finite increasing sequences (ti)ie{o, ....«} in 
[so, si] with Iq = sq and t„ = si. 

If 7 is a path [sq, si] — > X and 7' is a path [si, S2] ^ ^ such that 7(si) — 
7'(si), then we define their concatenation 77' to be the path [sq, S2] X that 
coincides with 7 on [sq, si] and with 7' on [si, S2]. 

The following result is a consequence of |(SJ). 

Lemma 3.5. For all paths 7 : [so,si] — > X and 7' : [si,S2] — > X smc/i that 
l{si) = 7'(si), we have 

= /(7) + Hi) ■ 

Assumption 3.6. We assume that for all t and x and y in X , A*{x,y) = 

sup^{/(7)}, where the supremum is taken over all paths 7 : [0,t] ^ X from x 
to y. 

This assumption can be shown to be satisfied in the case when A* = S**, the 
Lax-Oleinik semigroup of example 12.21 However, 1(7) will not necessarily be 
the integral of the Lagrangian along 7. 

For any path a : R_(- X, define 

J„(s, t) := A*{a{0), + /(«![,,*]) - A*{a{0),a{t)) . 

Clearly, Jq, is always non-positive and Ja{s,t) — Ja{s,u) + Ja{u,t) for all s, u, 
and i in R with s < u < t. Using this observation, we can prove the following 
lemma. 

Lemma 3.7. Let 7 be an e-almost~geodesic, for some e > 0. Let y Cz X and 
e' > 0. Then there exists an e' -almost-geodesic 7' starting at y and a constant 
A G M such that j'{t) — j{t + A) for t large enough. 

The following lemma is elementary. 

Lemma 3.8. Let u be a super-harmonic function. Then 

u = sup {^luiw) +w) ^ sup ifJ-uiw) + w) . 

Proposition 3.9. 1/7 : M+ X is an almost-geodesic, then K{-,'-f{t)) con- 
verges as t tends to infinity to some w € 

Proof. Let e > 0. By Lemma IS77I we lose no generality by assuming that 7 is 
an e-almost-geodesic starting at b. So, for all s and t in M_|_ with t>s, 

A*{b, 7W) < e + A*{b, 7(s)) + A*(7(s), 7W) 
<e + A*ib,^{t)) . 

Therefore 

-e<A*{b,^{s))+K{j{s),j{t)) . (9) 
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Since ^ is compact, it suffices to check that all convergent subnets of i4r(-, 7(t)) 
have the same limit and that this limit is in Let {j{td))deD and {'j{se))eeE 

be two subnets of (7(t))teE_|_ , such that the nets 

{K{-,j{td)))deD and {K{-,j{se)))eeE 

converge respectively to some w and w' in Applying with s = Se and 
t = tdi and taking the limit with respect to d, we obtain —e < A* {b,^{se)) + 
w(7(se)). Taking now the limit with respect to e, we get — e < H{w',w). Since 
e is arbitrary and H{w',w) is non-positive, we obtain H{w',w) = 0. From 
Lemma 13.81 we deduce that w > ^,yj(w') + w' — H{w',w) + w' = w' . By 
symmetry, we conclude that w — w' , and so H{'w, w) — 0. 

Let y G X and e > 0. By Lemma [3.71 there exists an e-almost-geodesic 7' 
starting at y and such that K{-,j'{t)) converges to w. For all s and t in IR+ 
with t > s, 

A^y,j'{s))+A*{j'{s),j'{t)) > A*{y,j'{t)) - e. 

Subtracting A*{b,y{t)), and letting t tend to infinity, we get 

A"iy, 7'(s)) + w(7'(s)) > wiy) - e. 

Since e is arbitrary, we have w{y) < sup^g^^ z) + 'w{z)) — {A^w){y). Since 
any element of ^ is superharmonic, we get the equality. □ 

Lemma 3.10. Let h : X Kmax be harmonic and let x € X . Then for each 
e > 0, there exists a path 7 starting at x that is e- almost-optimal with respect 
to h. 

Proof. Fix a sequence e„ in (0, 00) such that X^i^o ^" ^ ^/^- Since h is harmonic, 
we can construct a sequence (a:„)„gN in X starting at := x such that 

h{xn) < e„ + A^{xn,Xn+i) + h{xn+i) for all n e N. 

Adding these inequalities, we get 

n — 1 n— 1 

h[x) <^A^[xi,Xi+i) + h{xn) + for all n e N. (10) 

i=0 j=0 

By assumption 13.61 for each i € N, we can find a path 7^ : [i, z + 1] X with 
7i(i) = Xi and 7i(z + 1) = Xi+i such that 

/(7i) > 2:^+1) - e» (11) 

We define 7 : IR+ ^ X to be the concatenation of these paths, meaning that 
7(f) := 7[tj (t) for all t G M+. Here [ij denotes the largest integer no greater 
than t. To show that 7 is an e-almost-optimal path with respect to /i, consider 
any increasing sequence of times to, ■ ■ ■ ,tm with to '■— 0. Let n := \tm\ be the 
least integer no less than Adding the inequalities obtained from by 
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letting i run from to n — 1, and then using (|10|l . we get 

= /(7o) + --- + /(7n-i) 



n-1 



n-l 

>h{x)-h{xn)~2'^ej. (12) 

Using the fact that A""*™ (7(t,„), .x„) < h{j{t„i)) — h{xn), we get 
/(to, . . . ,tm;7) = /(to, . . . ,^m,f^;7) - A""*" (7(t„), a;„) 

> ^(7l[0,n]) + h{^n) - h{j{trn))- 

The e-almost-optimahty of 7 foUows on combining this with 1)12(1 . □ 

The following is a max-plus analogue of the Martin representation theorem 
in potential theory. A discrete time version was given in ;6 . 

Theorem 3.11. A function h : X ^ IRmax is harmonic if and only if it can be 
written 

h — sup iv{w) + w) , (13) 

where v is some upper semicontinuous function from ^™ to Rmax- Moreover, 
fih is the greatest v satisfying this equation. 

Proof. Let h be harmonic. By Lemma 13.81 

h — sup {iih{w) + w) > sup (/i/i(w) + w) . 

To show the opposite inequality, fix x € X such that h{x) > —00. Let e > 0. By 
Lemma 13.101 there exists an e-almost-optimal path 7 : — > X with respect to 
h starting at x. By Lemma 13.31 this path is an almost-geodesic and therefore, 
by Lemma [3.91 K{-,^{t)) converges to a point ^ of . Also, the e-almost- 
optimality of 7 implies that 

h{x)<e + A*{x,^{t)) + h{-i{t)) 

= e + K{x, 7(0) + A*{h, j{t)) + h{j{t)) . 

Letting t tend to infinity and using the fact that e is arbitrary, we see that 
h{^) — ^{^) + t^hiO- We have thus estabhshed that 

h{x) < sup {iih{w) + w) . 

Now let V : Mmax be any function satisfying and let ^ G 

Then h > + ^ So 

fJ-hiO = hmsup {A*{b,x) + h{x)) 
> lim sup {A* {b, x) + ^{x)) + 

= . 
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Therefore v < ^h- 



□ 



The remaining resuhs can be proved by adapting the arguments of the cor- 
responding results in 6^. 

Lemma 3.12. Let ^ e ./# be such that H{^,^) — and suppose that ^ can be 
written in the form ^ — sup^^^{iy(w) + w), where v : ^ Rmax is upper 
semicontinuous. Then, there exists w € ^ such that ^ = v{w) + w. 

Definition 3.13. Let y be a set of functions from X to M that is closed under 
maximisation and the addition of a constant. A function / from X to R is said 
to be normalised if f(b) ~ 0, and is said to be an extremal generator of V if 
/ = max(u, v) with u and v in V implies that / is equal to either u or v. 

Corollary 3.14. Every element ^ of ^ satisfying H{^,^) = is a normalised 
extremal generator of .y . 

Theorem 3.15. The normalised extremal generators of Jif are precisely the 
elements of . 

Proposition 3.16. Assume that ^ is first countable. For all ^ e ./#™, there 
exists an almost-geodesic converging to f . 



4 The Horoboundary of a Linear Quadratic Con- 
trol Problem 

We next study the eigenproblem ©, for the Lax-Oleinik semigroup (A*)t>o 
corresponding to the Lagrangian |0J. 

Observe first that A*0(O) > (j){Q), since the zero control makes the trajectory 
stay at the origin with a zero cost. Hence, any eigenvalue A of the semigroup 
must be nonnegative. 

The normalised semigroup Aj^ = — AT + A'^ corresponds to the reward: 

'^(|x(t)|2 + |i(t)P + A)dt, 

meaning that Aj^ixjy) is the maximal value of this reward over all curves x 
joining x to y. 

One may regard the optimisation of this reward as a physics problem in which 
a particle is subjected to a force away from the origin proportional in strength 
to its distance from the origin. Indeed, the Euler equation gives that optimal 
paths satisfy ic{t) = x(t) for all t S M. The general solution to this differential 
equation is x(i) = We* + Ze~*, where W and Z are constant vectors. If the 
position is x at time and y at time T, then 

W^y^^ and Z= 



Doing the integration, we find that the action of the optimal path is 

i\x\^ + \y\')coshT^2x-y 

^nhf ■ ^^^^ 
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Now we must maximise this expression over T > 0. We will first do the 
case when A = since the formulae are more manageable. If a; • y < 0, then 
the supremum occurs at T = +cxd, where the function above tends to the limit 
A*{x,y) = — — On the other hand, if x • y > 0, then the supremum is 
attained at finite T. Differentiating the expression above and setting it equal 
to zero, we obtain a linear equation in coshT, with solution 

coshT.^i^. 

2x ■ y 

From this, we can easily calculate sinhT as well. Substituting these values 
into H14|) and simplifying, we obtain 

A*{x,y) = -\x ~ y\\x + y\ iix-y>0. 

To calculate the boundary points we take the limit of A*{x,y) — A*{0,y) as 
\y\ tends to infinity, keeping n := y/\y\ constant. One expands the expression 
in powers of \y\ and discards those that tend to zero. The resulting boundary 
point is the function: 



hn{x) 



2{x ■ nf, if X • n > 0, 
otherwise. 



Figure n shows the "horospheres" of /i„ in dimension two, when n — (0, 1). 
These are the sets of the form {x g | hn{x) = const.}. 

Now we turn to the case when A > 0. To find the maximum over T > 
of l|14|l . we again differentiate it and set it equal to zero. This time we obtain a 
quadratic equation in coshT, with solutions 



AcoshT ^~x-y± y/ {x ■ y)^ + + X{\x\^ + \y\^) . 

It is of course the positive solution we want. We next calculate sinhT and 
T. Substituting them into (|14|l . we obtain an expression for A*{x,y), which, 
unfortunately, is now rather complicated. However the special case when a; = 
is not too bad: 



A*{Q,y) := ~\y\V^ ~ WogiVX + lf + \y\) + AlogA. 

To calculate the boundary points we again take the limit of A*{x, y) — A*{0, y) 
as lyl tends to infinity, keeping n := y/\y \ constant. The best way of calculating 
this limit is to treat separately the terms involving a logarithm from the rest. 
Again one should expand the expression as powers of \y\, discarding those that 
tend to zero. After much simplification, one arrives at the result 

M.-) = -A^+-"^^-Alog-^, 



where R :— ^ [x ■ n)^ + A — x ■ n. 

Figure El shows the horospheres of /i„ when n — (0, 1), again in dimension 
two. 



11 



Figure 1: Horospheres of the LQ model when A = 0. 
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